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Current-driven vibrational non-equilibrium induces vibrational sidebands in single-molecule tran- 
sistors which arise from tunneling processes accompanied by absorption of vibrational quanta. Unlike 
conventional sidebands, these absorption sidebands occur in a regime where the current is nominally 
Coulomb blockaded. Here, we develop a detailed and analytical theory of absorption sidebands, in- 
cluding current-voltage characteristics as well as shot noise. We discuss the relation of our predictions 
to recent experiments. 



PACS numbers: 73.23.Hk,73.63.-b,85.65.+h 



I. INTRODUCTION 

Single-molecule junctions and transistors have opened 
new vistas for quantum transport phenomena in nanos- 
tructures, which are associated with the coupling of 
the electronic degrees of freedom to collective modes 
of the molecule such as a molecular spin^'^ or molecu- 
lar vibrations.'^ One of the most fascinating possibilities 
afforded by this coupling is that the current flow can 
drive the collective modes far out of thermal equilibrium, 
which will in turn act back on the current. This makes 
single-molecule transistors an interesting testbed for out- 
of-equilibrium quantum transport. Indeed, several recent 
experiments^'^'^ on single-molecule transistors have pro- 
vided convincing evidence for out-of-equilibrium vibra- 
tions and theory has made numerous prcdictions^'^'^ for 
new effects arising from vibrational nonequilibrium. 

In the Coulomb blockade regime of transport through 
a single-molecule transistor, the coupling of the tunnel- 
ing electrons to molecular vibrations induces vibrational 
sidebands likiS, These sidebands occur for bias and gate 
voltages, where the Coulomb blockade is lifted, and arise 
from sequential tunneling processes which are accompa- 
nied by excitations of the molecular vibrations, cf. Fig. 
[IJa). If we consider tunneling processes which make the 
molecule switch between its neutral and singly-charged 
state, vibrational sidebands occur whenever EF + eV/2 = 
ed + nhio due to tunneling into the molecule and = 
Ep — eV/2 + nfvjj due to tunneling off the molecule. Here, 
Ep is the Fermi energy of the leads in the absence of 
an applied bias V , ed denotes the energy of the relevant 
molecular orbital, and to is the frequency of the vibra- 
tional mode. The positive integer n labels the different 
vibrational emission sidebands, which typically weaken 
with increasing n. 

These sidebands occur due to electron tunneling pro- 
cesses which are accompanied by emission, i.e. excita- 
tion, of vibrational quanta. In the presence of current- 
driven out-of-equilibrium vibrations, it is natural to 
ask whether there exist additional features in transport 
through single-molecule transistors which originate from 
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FIG. 1: (Color online) Schematic vibrational sidebands in 
dl/dV vs. bias and gate- voltage plane for (a) fast, (b) in- 
termediate, and (c) slow relaxation. Conventional emission 
sidebands in the sequential tunneling region (colored/shaded) 
are shown in white, absorption sidebands in red/gray. Hori- 
zontal dashed lines appear due to inelastic cotunneling, thick 
black lines due to sequential tunneling without change of the 
vibrational state. 



the absorption of vibrational quanta (similar to anti- 
Stokes lines in Raman spectroscopy). Indeed, it was 
recently observedii in numerical simulations of trans- 
port through single-molecule junctions that vibrational 
non-equilibrium can induce additional absorption side- 
bands within the region where the current is nominally 
Coulomb blockaded. A schematic illustration of these 
absorption sidebands is shown in Fig. [TJb) and (c) . 

It is the purpose of this work to present an analytical 
theory of these novel vibrational absorption sidebands. 
In Secini we discuss the basic mechanism which leads to 
the formation of absorption sidebands and identify sev- 
eral qualitatively different regimes. Our analytical the- 
ory, including both coductance and shot noise, is pre- 
sented in Sec. IIIII In Sec. llVi we discuss the relation of 
our results to recent experiments^!^ and argue that ob- 
servation of the vibrational absorption sidebands allows 
one to measure the vibrational relaxation time in single- 
molecule transistors. Some technical details have been 
relegated to an Appendix. 
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FIG. 2: Relevant transport processes in the vicinity of the 
first absorption sideband for huj < eV < 2huj. In addition to 
elastic cotunneling (not shown), the following processes con- 
tribute to current flow: (a) Inelastic cotunneling with emis- 
sion of one vibron; (b) population of the molecule by inelastic 
sequential tunneling with absorption of one vibron; depopula- 
tion by (c) elastic sequential tunneling or (d) inelastic sequen- 
tial tunneling with emission of at most two vibrons. Molecular 
states are symbolized by |n, g) where n refers to the charge 
state and q to the vibrational state. 



II. BASIC MECHANISMS 

The absorption sidebands occur for gate and bias volt- 
ages where the current is nominally Coulomb blockaded. 
In other words, the energy level of the molecular orbital 
lies outside the bias window. This is illustrated in Fig. [2] 
where we depict the processes relevant in the vicinity of 
the first absorption sideband for fiio < eV < 2huj. 

The basic process which allows current to flow through 
the molecule is cotunneling in which an electron tunnels 
between the leads with intermediate virtual occupation 
of the molecule. Cotunneling necessarily leaves the vibra- 
tional state unchanged for bias voltages eV < hoj (elastic 
cotunneling). Assuming temperatures T fku, vibrons 
can be excited once eV > fku [inelastic cotunneling, cf. 
Fig. [5] (a)]. In plots of the differential conductance of 
the single-molecule transistor, inelastic cotunneling due 
to excitation of vibrons leads to features at eV = nhiu 
with n = 1, 2 . . ., independent of the gate voltage^^ cf. 
Fig. [II 

If the vibrons excited by inelastic cotunneling decay 
only slowly, their presence opens a new transport channel 
which is responsible for the occurrence of vibrational ab- 
sorption sidebands. Indeed, absorption of these vibrons 
by an electron may energetically allow sequential tunnel- 
ing processes onto or off the molecule even though they 
are forbidden in the absence of vibron absorption [cf. Fig. 
[5] (b)] . For parameters in the vicinity of the coexistence 
of the neutral and the singly-charged molecule, the nth 
vibrational absorption sideband occurs when sequential 
tunneling becomes possible with absorption of nfuv. Set- 
ting Ep = 0, the locations of the absorption sidebands 
are eT^/2 -I- nfko = Cd for tunneling onto the molecule 



(with ed > eV/2) and -I- nhio = —eV/2 for tunneling 
off the molecule (with < —eV/^). In addition, the 
condition eV > fuo must be satisfied because vibrons can 
be absorbed only after they are first excited in inelastic 
cotunneling events. 

Rates for allowed sequential tunneling processes are 
lower order in perturbation theory in the molecule-lead 
hopping amplitude than electron cotunneling. Thus, 
these sequential tunneling rates Wsoq are generally much 
larger than the cotunneling rates W^cot- For this reason, 
we expect that (if energetically allowed), inelastic cotun- 
neling events are rapidly followed by a sequential tun- 
neling event involving vibron absorption. This scenario 
implies that three different regimes need to be distin- 
guished: (i) fast vibrational relaxation: 1/r ^ VFsoq 3> 
Wcot; (ii) intermediate relaxation: VFsoq ^ l/r 3> Wcot; 
(iii) slow vibrational relaxation: VFsoq ^ Wcot ^ 
There are no absorption sidebands for fast vibrational 
relaxation. Here, we focus on the other two regimes 
of intermediate and slow vibrational relaxation, both of 
which exhibit such additional sidebands. Interestingly, 
there is a qualitative difference between (ii) and (iii), cf. 
Fig. [ijb) and (c): For intermediate relaxation, emitted 
vibrons will decay before the next inelastic cotunneling 
event, even in the absence of a sequential tunneling pro- 
cess. This implies that a sequentially tunneling electron 
can absorb nhw only if the preceding inelastic cotunneling 
event excited the molecular vibrations by nhw. As a re- 
sult, the nth absorption sideband occurs only for voltages 
eV > nhw. In contrast, for slow vibrational relaxation, 
the vibronic state can be excited to levels beyond eV by 
repeated inelastic cotunneling. Thus, higher absorption 
sidebands become visible as soon as inelastic cotunneling 
becomes possible for eV > huj. 

This picture of transport in the regime of vibra- 
tional absorption sidebands has several additional conse- 
quences, (i) With increasing vibrational relaxation time, 
the current through the junction will increase since each 
inelastic cotunneling event will be followed by an increas- 
ing average number of sequentially tunneling electrons, 
(ii) The fact that inelastic cotunneling events are rapidly 
followed by sequentially tunneling electron(s) effectively 
implies electron bunching and thus super-Poissonian shot 
noise, (iii) This bunching becomes stronger for larger 
electron- vibron coupling A. Therefore, vibrational ab- 
sorption sidebands become more prominent with increas- 
ing A. 

Vibrational absorption sidebands are amenable to an 
analytical description for both weak and strong electron- 
vibron interaction, A <C 1 and X ^ li^ Results for 
these regimes give good qualitative insight also for the 
crossover region of intermediate coupling constants A ~ 
1 , where an analytical description becomes more cumber- 
some. 
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III. QUANTITATIVE THEORY 
A. Model and basic processes 

In view of the generic nature of the phenomenon under 
discussion, our considerations are based on the Anderson- 
Holstein mode &^^'^^i^^ for a single-molecule device in 
which transport is dominated by one gate-tunable spin- 
degenerate molecular orbital coupled to a single har- 
monic vibrational mode with frequency uj, as described 
by the Hamiltonian 

— ^dnd + U Ud^Udi + fuvb^b + Xhw{b^ + b)nd- (1) 

Here, U is the charging energy, the operator da annihi- 
lates an electron with spin projection a on the molecule, 
and Ud ~ rid] + ridi ~ X^o- '^ct'^c denotes the corre- 
sponding occupation-number operator, b annihilates the 
vibrons. Tunneling between the molecule and (Fermi- 
liquid) source and drain electrodes with Fermi energies 
^iL,R = ±eT//2 is assumed to be weak with F ^ T. (F 
is the tunnel broadening of the molecular orbital in the 
absence of the vibrational mode and T the temperature.) 

Explicit results for the tunneling rates can be obtained 
by expanding the T-matrix in the molecule-lead tunnel- 
ing and applying Fermi's golden rule. To be specific, we 
focus on absorption sidebands in the Coulomb-blockade 
diamond where Ud — Q- Cotunneling processes between 
leads a and b involving a change of the vibrational state 
from q to q' occur with a rate (assuming U — > oo) 
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Here, s'^"* = 2 denotes a spin factor and /a(e) the Fermi 
distribution of lead a. Analytical results for the cotun- 
neling rates are obtained using the integrals 
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Here, the poles have been regularized in a standard 
fashioniiiiii We used the notation Ef- = ^ ± — ej] , 

the digamma function ^l^{z), and the Bose distribution 
nsiE). The Franck-Condon matrix elements Mq^q^ — 

(qi je'''^''^''^' 152) between vibrational states \q) have the 
analytical expression 
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in terms of generalized Laguerre polynomials L'^^^x) 
where q = min(qi,g2), Q = max(gi,g2)- Vibration- 
assisted sequential tunneling onto the molecule involving 
lead a and changing the vibrational state from q to q' has 
the rate 
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with the spin factors s^"*" = 1 and s°^^ = 2. An anal- 
ogous expression holds for W^j^,.^. Below, we will drop 
the index "ab" on the cotunneling rates from source to 
drain (LR). We also abbreviate the total cotuncling rate 
as %'°'tot = Ea,fc^gg';a6- Sequential tunneling rates 
without lead index refer to the sum of the rates across 
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B. Intermediate vibrational relaxation 

1 . Current 

For weak electron- vibron coupling A 1, crucial sim- 
plifications occur due to the Franck-Condon matrix ele- 
ments favoring processes with a minimal change in the 
vibron occupation. This allows us to restrict attention to 
the occupation probabilities of a limited number of 
states |n, q) of the molecule. Here, n denotes the charge 
state and q the vibrational state. For example, consider 
the system with A ^ 1 after an electron has sequen- 
tially tunneled into the molecular state |1,0) by vibron 
absorption. Tunneling out to the drain electrode can now 
be accompanied by emission of one or several vibrons [cf. 
Fig. [DJd)]. However, to lowest order in A, we only need 
to include the process with the final state | 0, ) in which 
no vibron is excited [cf. Fig.[2Jc)]. 

We first focus on the nth vibrational absorption side- 
band, for eV > nhcu and intermediate relaxation. The 
states |0,0), \0,j < n), and |1,0) need to be included to 
leading order. Solving the corresponding rate equations, 
we obtain for the stationary current (see appendix |^ for 
details) 
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The first term on the right-hand side accounts for con- 
tributions of elastic and inelastic cotunneling. The sec- 
ond term is responsible for the nth absorption sideband. 
Here, the inelastic cotunneling rate Wq°^ reflects that the 
corresponding process opens the absorption-assisted se- 
quential tunneling channel by which additional electrons 
can pass through the molecule. The first fraction arises 
from the competition between sequential tunneling and 
vibrational relaxation. It is the rate W^q ^ which turns 
on at the nth absorption sideband and thus leads to a 
peak in the differential conductance dl/dV vs. bias and 
gate voltage when eV /2 + nhuj = td- The second fraction 
describes that an additional electron passes through the 
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molecule only when the electron on the molecule sequen- 
tially tunnels out to the right lead (drain). 

For strong electron- vibron coupling A » 1, the Franck- 
Condon matrices favor transitions where the final state 
has the largest accessible vibron occupation. Unlike for 
weak coupling, this allows for longer sequences of more 
than two sequential tunneling events. This is because 
tunneling out from the molecule will now preferredly lead 
to final states with vibronic excitations [e.g., process (d) 
in Fig. [2] will be dominant with respect to process (c)]. 
Hence, the restriction to only a small set of states is no 
longer feasible. Instead, the system can acquire a large 
number of vibrons. The sequence of sequential tunneling 
events is ultimately cut off by subdominant sequential 
processes or vibrational relaxation. In appendix |^ this 
is illustrated in detail for the first vibrational absorption 
sideband huj < eV < 2huj where we find 
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Here, the first term describes the elastic cotunneling and 
the second term the inelastic cotunneling which initi- 
ates the series of sequential events described by the third 
term. The expression Eq. ([5]) makes it again explicit that 
the first vibrational sideband emerges for intermediate 
relaxation, i.e., once the sequential tunneling rate W\ 
becomes faster than the vibrational relaxation 1 /r 
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FIG. 3: The quantity F = l-i- hcq/Icot as defined in the text 
vs. the relaxation time for (a) weak (A = 0.2) and (b) strong 
(A = 4) electron- vibron coupling. {V = 1.7 huj; ea = 1.6 huj) 



the nth absorption sideband with eV > nHuj, we obtain 
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Due to the first fraction, the Fano factor for weak cou- 
pling will take a value only slightly above one. 

Much larger Fano factors are to be expected for strong 
electron-vibron interaction where cotunneling events can 
be followed by a long sequence of vibron assisted sequen- 
tial tunneling events. An approximate measure of the 
Fano factor is given by F = l-l-/scq//cot, where Iscq (/cot) 
denotes contribution to the current due to sequential (co- 
tunneling) processes. This is an instructive measure for 
the average number of sequential tunneling events follow- 
ing a cotunneling process. We plot numerical results for 
F as a function of \ogiQ{TT/h) for both weak couphng 
[Fig.IHa)] and strong coupling [Fig.[3ljb)]. In both cases, 
F features a step-like increase at the crossover between 
fast and intermediate relaxation but the step height is 
much greater for strong electron-vibron coupling. (Note 
the different scales on the vertical axes in Fig. [3l) 



2. Fano factor 



Slow vibrational relaxation 



An interesting feature of Eqs. ([7]) and ^ is that the 
current is partially vibron assisted in the sense that it 
increases with increasing vibrational relaxation time r. 
This is because the vibron-assisted sequential tunneling 
processes which may follow inelastic cotunneling event, 
transfer additional electrons. Clearly, the number of ad- 
ditional electrons passing through the molecule increases 
when going from fast to intermediate vibrational relax- 
ation. E.g., for A <C 1 the average number of electrons 
passing per inelastic cotunneling event is 1 for fast relax- 
ation, but 3/2 (symmetric junction) or 2 (strongly asym- 
metric STM configuration) for intermediate relaxation. 

Since sequential tunneling proceeds on a much faster 
timescale, these considerations imply electron bunching 
and super-Poissonian shot noise. For A <C 1, the cor- 
responding Fano factor F can be obtained from F — 
{Ii + 2/2)/(/i + ^2), where h {I2) denotes the contribu- 
tion to the current due to single-electron (two-electron) 
processes. For intermediate vibrational relaxation near 



For voltages eV < nhuj, the nth vibrational absorption 
sideband becomes visible only for slow vibrational relax- 
ation, since the molecular vibrations must be excited to 
the nth level by repeated inelastic cotunneling processes. 
To make this point quantitative, we consider the second 
absorption sideband for voltages hoj < eV < 2Huj and 
weak electron-vibron coupling A 1. (Analytical re- 
sults for strong coupling become more cumbersome and 
will not be discussed here.) We find (see appendix [X\ for 
details) 
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where I^ot = ePl^W§§^ + eW§'^^ Pl^[l + W^l^y D + W^f / D] 
collects cotunneling processes without subsequent se- 
quential tunneling, D = Wf^^ + PFfo'tot + 1/^, and 
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The product of two cotunneling rates appearing in the 
second term of Eq. pil)l reflects the need for two sub- 
sequent inelastic cotunneling events to excite the vibra- 
tional level q~2. Due to the denominator D the second 
absorption sideband appears for these voltages once vi- 
brational relaxation becomes slow compared to inelastic 
cotunneling. 



IV. CONCLUSIONS 



the Coulomb blockade region, where several "molecular 
orbitals" are located within the bias window. 

Finally, we remark that the magnitude of the vibra- 
tional absorption sidebands may also allow access to the 
"effective temperature" of the out-of-equilibrium molecu- 
lar vibrations which is an important characteristic of the 
single-molecule junction but has so far remained outside 
experimental reach. 



Vibrational absorption sidebands exhibit a rich phe- 
nomenology due to the presence of three relevant time 
scales derived from the rates for cotunneling, sequen- 
tial tunneling, and vibrational relaxation. This leads to 
the existence of different regimes, and as a consequence 
the observation of vibrational absorption sidebands can 
be employed to measure the vibrational relaxation time. 
This is particularly promising in scanning tunneling mi- 
croscope measurements where the current level can be 
tuned by changing the tip-to-molecule distance. In this 
case, the vibrational absorption sidebands should emerge 
only once the current exceeds a certain threshold current, 
and one could extract the vibrational relaxation time r 
in two different ways: 

(i) While vibrational absorption sidebands emerge al- 
ready when 1/t « Wsoq, one suggestion is to focus on 
higher-order [n > 1) absorption sidebands for hw < 
eV < 2huj. These sidebands appear only once the vi- 
brational relaxation rate becomes comparable with in- 
elastic cotunneling rates. Assuming electron-vibron cou- 
pling constants of order unity, the current is given by 
the inelastic cotunneling rate (multiplied by e) up to a 
numerical factor of order unity. If these sidebands ap- 
pear at a threshold current of order Ic, the vibrational 
relaxation rate can be estimated from l/r ~ Ic/e- For 
electron-vibron coupling constants A which are not too 
close to unity, a more accurate quantitative comparison 
with the theory developed here should provide a good 
estimate for l/r. 

(ii) While suggestion (i) requires slow vibrational re- 
laxation, an alternative approach which applies also to 
intermediate vibrational relaxation consists of a separate 
measurement of the sequential tunneling rate Wseq- This 
is readily done by measuring the current at the same bias, 
but at a gate voltage such that the single-molecule junc- 
tion is outside the Coulomb-blockaded regime, where the 
current is dominated by sequential tunneling. 

Indeed, recent experiments^ have extracted vibra- 
tional relaxation times for suspended single-wall carbon 
nanotubes in a related way and found them to be of the 
order of 10ns. However, it is important to understand 
that the regime discussed in the present work differs from 
that investigated in these experiments in an essential way. 
The vibrational absorption sidebands discussed here oc- 
cur for gate and bias voltages where the current is nomi- 
nally Coulomb blockaded. By contrast, the experiments 
of Refs. iJ andd are performed at high bias, well outside 
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APPENDIX A: CALCULATION OF THE 
CURRENTS FROM REDUCED RATE 
EQUATIONS 

Here, we sketch the derivations for the current expres- 
sions Eqs. (O, (O, and pUj) . valid in different regimes of 
relaxation strength and electron-vibron coupling. The 
general procedure is as follows: For each particular 
regime, we identify the relevant molecular states and 
transitions between them. The corresponding rates are 
calculated from a T-matrix expansion as discussed in sec- 
tion [111X1 Many rates are found to be essentially zero, 
either because the process is not allowed at small temper- 
atures (T ^ eVyhio), or because the vibronic transition 
is suppressed due to the Franck-Condon matrix element. 
Eliminating these processes reduces the system of rate 
equations determining the stationary occupation proba- 
bilities according to 

dP" r n 

= — - = 'V P"] " - P"W^"? 

[ q' "g'q q '^^ qq' 

n' ,q' 

-^Pq-Pq^T^q^Pq']- (^1) 

For simplicity, we have used the notation W™i for co- 
tunneling rates. The last term in Eq. (|Aip models re- 
laxation towards the equilibrium vibron distribution P^'^ 
with time constant r. The stationary current is then 
given by the sum of the contributions from sequential 
tunneling and cotunneling, 

n,q,q' 

+ eJ2p- [W^l^j,, - W-^.^„^j] . (A2) 

n,q,q' 



FIG. 4: Relevant states and processes for intermediate re- 
laxation and weak electron-vibron coupling in the vicinity of 
the nth absorption sideband. Single solid lines represent co- 
tunnehng transitions, double lines sequential tunneling and 
dashed lines vibrational relaxation. 



FIG. 5: Dynamics leading to the first absorption sideband (in 
the bias range hco < eV < 2hjjj) for intermediate relaxation 
and strong electron-vibron coupling. 



1. Intermediate relaxation, weak electron-vibron 
coupling 

We first consider the nth absorption sideband for A ^ 
1 and intermediate relaxation. The states and processes 
relevant in this regime are depicted in Fig. |4l The left 
and upper right branches of the scheme (originating from 
the ground state) represent elastic and inelastic cotun- 
neling channels. These compete with the series of pro- 
cesses (lower right branch) which is responsible for the 
formation of the nth sideband. The reduced system of 
rate equations which describes this interesting part of the 
scheme is 
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Here we have used Eq. (IA1[) with Pq^ « 1, valid for 
T hio. ^From this we obtain the probabilities 
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Insertion of Eq. (IA5p into the equation for the current, 
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with P(? 



1 - Wcot/(M^soq + 1/t) « 1 yields Eq. 0. In- 
tuitivcly, the last approximation can also be understood 
from the fact that the system will predominantly occupy 
the ground state. While inelastic cotunneling makes the 
system leave the ground state, such excursions only hap- 
pen on a very short time scale, limited by the quick re- 
turn to the ground state via sequential tunneling or vi- 
brational relaxation. 



2. Intermediate relaxation, strong electron-vibron 
coupling 

We now turn to the limit A ^ 1 in the intermedi- 
ate relaxation regime. We focus on the first absorption 
sideband in the bias range hiu < eV < 2huj. The cor- 
responding dynamics is sketched in Fig. [51 Unlike in 
the A ^ 1 case, here, sequential tunneling out of the 
molecular orbital predominantly leads to highly excited 
vibronic final states. From there, the molecule can be re- 
populatcd by inelastic sequential tunneling. As a result, 
the initiating inelastic cotunneling event may be followed 
by a chain of several sequential processes. The reduced 
system, depicted in Fig. [5l includes not only processes 
of lowest order in 1/A but also higher order transitions 

1, q ) I 0, ) which serve as limiting processes for the 
chain of sequential transitions. By contrast, transitions 

1, q ) ! 0, g' ) with < g' < g -|- 2 are neglected since 
they produce only a higher-order correction to the se- 
quential current. The stationary current is 
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Here, the probabilities P/ are determined from the rate 
equations corresponding to Fig. [5l The resulting system 
of linear equations in this case is infinite. However, its 
simple structure allows for a recursive solution of in 
terms of P^+i, which can in turn be expressed in terms 
of Pl_i and so on. For strong electron-vibron coupling, 
we have W/°^_2.^ ^ ^lo-R^ ^^^^ second term in 



the sum in Eq. (|A7|) may be neglected. Approximating 
Po° « 1, we arrive at Eq. ([5]). 



3. Slow relaxation, weak electron-vibron coupling 

We finally turn to the regime of slow relaxation and 
weak electron-vibron coupling. We focus on the second 
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FIG. 6: Relevant states and transitions at the second absorp- 
tion sideband (in the bias range hijj < eV < 2hu!) for slow 
relaxation and weak electron-vibron coupling. 



absorption sideband in the bias range Hlj < eV < 2Huj. 
As shown in Fig. [SI at least two successive cotunneling 
transitions are required to open the sequential channel 
which generates the absorption sideband. The sketch is 
easily translated into a system of rate equations. Solving 
for the occupation probabilities and using Pq + P° « 1 
(up to corrections of order Wcot/W^soq) leads to Eq. (|10p . 



^ H.B. Heersche et al, Phys. Rev. Lett. 96, 206801 (2006). 

^ M.H. Jo et ai, Nano Lett. 6, 2014 (2006). 

^ For a recent review, see M. Galperin, M.A. Ratner, and A. 

Nitzan, J. Phys.: Condens. Matter 19, 103201 (2007). 
* B.J. LeRoy et al, Nature 432, 371 (2004). 
5 B.J. LeRoy et al., Phys. Rev. B 72, 075 413 (2005). 
^ D. Djukic, thesis (Leiden, 2006), |http:77hdrhandle[ 

net/1887/4927. 

A. Mitra, I. Aleiner, and A.J. Millis, Phys. Rev. B 69, 

245302 (2004). 
8 L.Y. Gorelik et al, Phys. Rev. Lett. 80, 4526 (1998). 
^ J. Koch and F. von Oppen, Phys. Rev. Lett. 94, 206804 

(2005). 

^° H. Park et al, Nature (London) 407, 57 (2000). 

J. Koch, F. von Oppen, and A.V. Andreev, Phys. Rev. B 
74, 205438 (2006). 



V. Aji, J.E. Moore, CM. Varma, |cond-mat/0302222| 
(2003). 

For A S> 1, it may happen that Wseq ^ Wcot is no longer 
valid, see Ref. [ll|. We do not consider this situation. 
" L.L Glazman and R.L Shekhter, Sov. Phys. JETP 67, 163 
(1988). 

N.S. Wingreen, K.W. Jacobsen, and J.W. Wilkins, Phys. 
Rev. Lett. 61, 1396 (1988). 

S. Braig and K. Flensberg, Phys. Rev. B 68, 205324 (2003). 
" D.V. Averin, Physica B 194-196, 979 (1994); M. Turek 
and K.A. Matveev, Phys. Rev. B 65, 115332 (2002). 
Strictly speaking, this statement holds true in the relevant 
range of transistions between low-lying vibrational states 
below the Franck-Condon parabola. 



